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Abstract. The patterns arising from the differential swelling of gels are investigated experimentally and 
theoretically as a model for the differential growth of living tissues. Two geometries are considered: a thin 
strip of soft gel clamped to a stiff gel, and a thin corona of soft gel clamped to a disk of stiff gel. When 
the structure is immersed in water, the soft gel swells and bends out of plane leading to a wavy periodic 
pattern which wavelength is measured. The linear stability of the flat state is studied in the framework 
of linear elasticity using the equations for thin plates. The flat state is shown to become unstable to 
oscillations above a critical swelling rate and the computed wavelengths are in quantitative agreement 
with the experiment. 

PACS. 46.32.4-x Static buckling and instability - 61.41. +e Polymers, elastomers, and plastics - 87. 18. La 
Morphogenesis 
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O ' 1 Introduction hand they were shown to be important in phyllotaxis (the 

^ ■ arrangement of leaves in plants) ^G] > m the wrinkling of 

Living organisms are full of fascinating complex patterns, leaves 00], in the selection of cell sizes jS], as well as in the 

One might wonder about the physical mechanisms at stake development of embryos On the other hand, a theo- 

as well as their relevance. Although a tissue is obviously retical framework was introduced to study the instabilities 

an elastic solid, the role of mechanical stresses in mor- occuring in the growth of elastic bodies . 
phogenesis was not recognised until recently. On the one 

At first sight, it is difficult to find physical systems al- 
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changes when submitted to external stimuli such as vari- We prepared our gels as in . A mixture of acry- 

ations in temperature, pH, osmotic pressure [HE], electric lamide (AA) and N,N'-methylenebisacrylamidc (BISAA) 

field ^U] or light Strictly speaking, such a gel does is dissolved with sodium acrylate (SA) in distilled water, 

not grow but swells by absorbing water while its elastic The polymerization is initiated by ammonium persulfatc 

modulus decreases. A number of studies (see ^1 f° r a (PA) and is catalysed with N,N,N',N'-tetramethylenediaminc 

review) have been devoted to the instabilities of swelling (TEMED) (0.3% in volume). The composition of the gels 

or deswelling gels, such as the folding of the surface of used mostly is given in Tab. ^ I n these conditions, gela- 

swollen gels clamped to hard substrates and the subse- tion generally occurs within one minute after the addition 

quent formation of a network of cusp lines 13,14,15,16 of the catalyst. 



HZj 



The characteristics of the gel can be tuned by varying 

In this article we are concerned with the instabilities 



the concentrations of the components. The more concen- 
occuring in the differential swelling of gels. Our motivation trated (and; for & same concen t ra tion, the more concen- 
is to design physical counterparts to growing tissues. We trated in BISAA) the so i utiorij the differ the gel. Likewise, 
use thin gel plates made by assembling two gels with differ- thc swelling rate can be incre ased by adding sodium acry- 
ent elastic and swelling properties and we investigate the late For the purpose of the experiment we have prepared 
resulting patterns experimentally and theoretically First, twQ digtinct types of gel . (T) a goft and swelling gel; (II ) a 
in section El we describe the experimental procedure and gtiff and nonswelling geL The elastic and swe m ng proper- 
give our first observations. In section we introduce the tieg of thege gelg haye been measured and are reported m 
theoretical framework and the linear stability analysis. In ^ 
section^] we compare quantitatively the experimental and 

theoretical results. We eventually give a discussion and j jj 
some perspectives in section |SJ 

2 The experiment 

Our goal was to design an experiment to study the differ- 
ential swelling of gels. We chose to assemble two thin flat 

gels having different elastic and swelling properties. Poly- Table 1 ■ Composition and properties of the gels used in the ex- 

, n ., ,i ii i j periment. Concentrations are given in mmol.L -1 . The swelling 
ry lamide gels were suitable as they swell when immersed 

rate corresponds to the increase in linear dimensions. E is the 



[AA+BISAA] 


720 


2880 


BISAA:AA ratio 


1:37.5 


1:19 


[SA] 


46-183 





Swelling rate 


50%-80% 


6% 


E (Pa) 


5.0 ■ 10 3 


3.2 ■ 10 5 



ac 



in water whereas their swelling and elastic properties can 
be tuned independently. 



elastic modulus of the gel. 
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Fig. 1. The strip geometry. Schematic, top and side pictures. 
A strip of soft swelling gel (I) is chemically clamped to a strip 
of stiff nonswelling gel (II). When immersed in water, part (I) 
swells while part (II) does not. An off-plane instability appears 



The gelation process is performed in a thin cell com- . * "k Soft swelling gel (I) 

posed of two glass plates separated by rubber spacers of 
constant thickness h (1 to 5 mm). The cell is set vertically. 
The spacers are used as masks shaped according to two 
different geometries: 

— The so-called strip geometry consists of two thin strips 
of gel of respective compositions I and II, clamped by 
their edges (see Fig.^l. To obtain this geometry, solu- 
tion II was first poured into the cell and left for gela- 
tion. Then solution I was added to form a second layer 

on (I) in the x direction if the width I is large enough (there are 

of width I (typically 1 cm) and length 20 cm. In the , , , jn . x T , , , . 

w x- j i o regularly spaced cusps otherwise) . Its wavelength A increases 

process of gelation of the second layer, the two layers monotonically with the width , as seen in pictures . 
become chemically clamped to each other. 

- The corona geometry is the axisymmetric counterpart ments with different swelling ratios and found similar re- 
of the strip geometry. A disk (radius in the range suits. For a small width I no buckling occurred but we ob- 
2-5 cm) of gel II is clamped to a corona (inner radius served instead a fine pattern with regularly spaced cusps 
r t , outer radius r a in the range 2-5 cm as well) of gel similar to those studied in [T2| . In the corona geometry, the 
I (see Fig. HJ. We used circular masks of various radii pattern is periodic and is characterised by a wavenumber 
to obtain these shapes: first the disk is made, then the (defined as the number of complete wavelengths) which 
cell is opened, the mask is replaced and the corona is increases with the aspect ratio r a /ri. A more quantitative 
moulded after closing back the cell. description is delayed to section g| 

The obtained structure is removed from the cell and 
immersed in water for several hours until a stationary o ti i . . ■ 

J 3 Theoretical setting 

state is reached and the swelling process is complete. In 

both cases, the sample undergoes a mechanical buckling We now study the patterns using a linear stability anal- 



instability such that the soft gel is no longer flat (except ysis of the equations of elasticity. We consider only the 
for narrow strips - see below) . This results in a wavy pat- soft swelling part (I) of the structure, the other part (II) 



tern (Figs.[IH2J with a well-defined wavelength. 



being considered as static as its swelling rate is small and 



In the strip geometry, the wavelength increases with it is very stiff. Since the gel sample is thin and flat be- 
the width of the strip I. We also performed some experi- fore the instability, we can use the Foppl - von Karman 
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Fig. 2. The corona geometry. A corona (inner radius r% and 
outer radius r ) of soft swelling gel (I) is clamped to a disk of 
stiff nonswelling gel (II). When swelling, the corona becomes 
unstable and goes off the inital plane. The wavenumber (eight 
on this picture) increases monotonically with the aspect ratio 
n/r - 



I -a 2 
E 



(Uyy + <JU XX ) , 



(3) 
(4) 



where a = 1/2 is the Poisson ratio of the gel (these gels are 
almost incompressible). The Foppl - von Karman equa- 
tions for equilibrium |19| read 



dx a 
da af 3 



= 0, 



(5) 
(6) 



for out of plane bending and in-plane stretching respec 
tively. The bending stiffness is 

Eh 3 



where D = 



(7) 



12(1 -o- 2 )' 

The corresponding (linearized) boundary conditions are 
given in the next two subsections. 



equations for thin elastic plates {H^. Material points are 
parametrized by their initial planar cartesian coordinates 
x and y. A deformation is defined by the displacement field 
(u x (x,y),u y (x,y),((x,y)); u x and u y are the in-plane dis- 
placements along the x and y axes respectively whereas 
£ is the transverse (off-plane) displacement. We use the 
framework of linear elasticity which will prove sufficient 
for the interpretation of the results (this restriction is dis- 
cussed in the conclusion), so that the in-plane stress ten- 
sor a a f3 (a, (3 = x, y) depends linearly on the deformation 
tensor u a p |19j : 



U aj 3 



_ 1 / du a dup d( d( 
2 \dx/3 dx a dx a dxp 

E 

\U XX ~r &W>yy) i 



1-cr 



(1) 
(2) 



3.1 The strip geometry 

We now specialise to the strip geometry. The reader should 
refer to Fig. ^for the notations. From this point onwards 
we use I as unit of length and D/(hl 2 ) as unit of stress. 
The constraint imposed by II results in the compression 
of the swollen strip in the x direction. We study the linear 
stability of the flat solution such that £ = 0. Using the 
translational invariance of the system we look for in-plane 
deformations in the form u x = —kx (k > 0), u y = f(y). 
The boundary conditions u y (0, x) = at the clamped edge 
and a xx (x, 1) = at the free edge, along with ©, impose 
u y = 0. Then the equilibrium equations reduce to 



A 2 ( - = 0. 
dx 2 



(8) 
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We look for periodic solutions in the form 

C(x,y) = ^(y)cosqx. (9) 
Using ©, appears to be the solution of a fourth-order Fig 3 Mogt ungtable mode for an infinite flat gtrip of gd 
linear differential equation with constant coefficients. Its i ongitu dinally compressed on one edge, as given by the solution 
solutions can be written in the form to ©. 

£(y) = Ae~xjp(imy)+B exp(— imy)+C exp(ny)+D enp(— ny), 

(10) 3.2 The corona geometry 

with 



(to 2 + q 2 ) 

The boundary equations on the clamped (y = 0) and 



k — V _LJ and n 2 — m 2 = 2q 2 (11) We now P rocee d to the same analysis for the corona ge- 
q 2 

ometry. The swollen corona (n < r < r a ) is supposed to 
be radially pulled at the inner edge so that its inner ra- 

free (y = 1) edges, read [T9*|: 

dius shrinks to fit the radius of the unswollen disk. Note 

£(o) = o, e"(i)-^ 2 e(i) = o, 

that Ti and r Q are the dimensions of the corona if it had 

£'(0)=0, -^{l) + {2-a)q 2 £ > '{\)=Q, 

swollen unconstrained. Although neither of these radii can 

These conditions can be viewed as a system of 4 linear 

be measured directly, their ratio Ti/r is the same as be- 

equations with the four unknowns A, B, C, D. A non-zero 

fore swelling. From this point onwards we use the outer 

solution £ exists if and only if the determinant is zero, 

radius r Q as unit of length and D/(hr 2 ) as unit of stress. 

which occurs for a certain n(q) (we recall that m can be 

In these reduced units, is the aspect ratio of the corona 

expressed as a function of to and q from Eq. Ijllll ). 

before the swelling. 

We now look for the most unstable wavenumber q c , i. 

. . . We first need the stress field in the flat corona as re- 

e. the one for which a non-zero solution exists with the 

. suiting from the tension at the inner radius. We look for 
lower value of the length increase k — k c . Ihus we lm- 

.... , , . . , , . . , . a solution of |JBJ in the form: 

phcitly assume that k is near its threshold value k c . lhis 

is another limitation which is discussed in the conclusion. ^ 

u r = ar H — , ug — 0. (13) 

The quantity k = (n(q) 2 + q 2 ) j q 2 is then minimized with 

respect to q. Eventually the wavelength reads (in dimen- f and are the standard polar coor dinates. The boundary 

sional units) conditions are u r (r = r,) = —/3 (displacement to fit the 



A 3.256 1. (12) unswollen gel) at the inner edge and a rr (r = 1) = (stress 



2tt 

A 3-dimensional representation of the corresponding solu- free) at the outer edge. The problem has now one single 
tion y) of JBJ is shown in Fig. degree of freedom /3 . We obtain a stress tensor of the 
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form: 

a rr = -a f 1 - , a ee =-a(l+ . (14) 

The number a being proportional to (3 is chosen as the 
swelling control parameter. The balance of moments JSJ 
reduces to 




^C + ^C + ^^ + -^ + ^j=o, (15) 

.,, . . <9 2 C 1 <9 2 C 1 Fie. 4. Most unstable mode for a flat corona of gel with 

with AC, = + ^r—; + - — . (16) 
or z r z o9 z r or 

radial tension at the inner radius. The aspect ratio (inner ra- 

At the clamped edge (r = fj), the boundary conditions 

dius/outer radius) is 0.74. The selected wavenumber is m — 8. 

read 

e = o, 

-g- = 0, (17) the conditions at the outer boundary (|18|) exists only for 



whereas at the free edge (r = 1), 

d I fd 2 C d 3 C 



a certain a(m,ri), which we compute numerically. Using 
the same argument as in the strip geometry, we choose 



dr r 3 \d9 2 drdO 2 J m(ri) for which a(m,ri) is minimum. Thus we obtain the 

39 2 dr 



1 fd 2 ( dC 

+ (& — 1)— 2 ( ~7T7y2 + r JT ) = 0- (18) wavenumber to as a function of the aspect ratio r^. Fig. 



r 

TTT , , . ... . . — _ . . r El shows the solution £(r,9) of (11511 with the aspect ratio 

We look for periodic solutions to (|15|) m the form 

Ti/r = 0.74 for which m = 8. 

C(r,0) =£(r)cosm0, (19) 

We also tested the r, — » 1 limit of our numerical cal- 

where the wavenumber m is an integer. We find a fourth- 

culation. In this limit, the corona geometry should reduce 

order linear equation in £(v). 

to the strip geometry, as the curvature of the interface 

We first compute a basis of the 2-dimensional space 

between the two gels vanishes. The wavelength reads 

formed by the solutions satisfying the conditions at the 
inner boundary i|17[) . To proceed, we solve the differential 

equation with initial conditions: \ n \ _|_ r 



l to 1 — r% 

(0,0,1,0) 



(20) 



(0,0,0,1) 

using a Runge-Kutta algorithm, so that we find the basis For j-j = 0.9999, we obtained indeed X/l = 3.255 in agree- 
£i, £2- A non-zero linear combination of £1 and £2 verifying ment with 1)12(1 . 
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Fig. 6. Wavenumber m of the instability as a function of the 
aspect ratio (inner radius / outer radius r ) before swelling. 
Crosses: experimental results. Line: theoretical result. 



Fig. 5. Instability wavelength A as a function of the width 

I of the swollen strip. The measurements are normalised by 4.2 The corona geometry 

the thickness h of the gel. Off-plane instability is observed for 
I > lo ~ 2h. For I < l c , cusped patterns appear on the surface 
of the swollen gel. The line corresponds to the theoretical result 
A = 3.256 1 (Eq.Ql- valid for I > h). 



4 Results 

4.1 The strip geometry 

The two geometrical parameters are the width I of the 
swollen strip, and the thickness h of the gel before swelling. 
We plotted our observations on Fig. [5] along with the 
analytical result i|12|l . The gel goes off the plane for I > 
l c ~ 2h. The wavelength of the instability increases with 
I, and is approximately linear with I in the limit of small 
thickness h <c I, in agreement with the theoretical result. 
For I < l c , instability patterns are observed on the surface 
of the swollen gel. They consist in regularly spaced cusps 
similar to those observed in 



In the corona geometry experiment, the thickness h was 
hxed to 1 mm, and the inner and outer radii were varied 
from 20 mm to 50 mm. Thus the theory of thin plates re- 
mains valid, which allows a direct comparison between ex- 
perimental and analytical results (Fig. |SJ). The wavenum- 
ber m increases with the aspect ratio Tijr , in quantitative 
agreement with the theoretical predictions. 

5 Conclusion 

To summarise, we showed that the swelling of thin soft gel 
plates clamped to a stiff gel leads to a buckling instabil- 
ity. A linear stability analysis yields a prediction for the 
pattern wavelengths in quantitative agreement with the 
experiment. These wavelengths are mainly determined by 
the in-plane geometry of the thin gel. Our analysis is re- 
stricted by two main limitations. 

On the one hand, we used the most unstable modes 
to predict the wavelengths at the instability threshold 
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whereas the experimental conditions are far above this References 
threshold. However standard weakly nonlinear analysis 
generally yields almost the same wavelengths [201 • On 
the other hand, we used the framework of linear elastic- 
ity while the swelling rate is large and the deformations 
are finite. Moreover the flat base state is anisotropic: for 
instance, in the strip geometry, the gel has swollen differ- 
ently along the two in-plane direction, so that the elas- 
tic modulus is not strictly the same in these two direc- 
tions. The small discrepancies between the experimental 
and theoretical wavelengths might be ascribed to either of 
these two limitations. 



A suggestion for future research stems from the sec- 
ond limitation. One might use the framework developed 
in [7j to build a more precise theory for the present exper- 
iment. On the experimental side, we have shown how to 
design physical counterparts to growing tissues. Other ge- 
ometries are currently investigated and will be the subject 
of forthcoming publications. 
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